library(dplyr) ; library(ggplot2) ; library(ggthemes) ; library(tidyr)
## Functions 
kable_cust = function(df_aux,ndec=5,cap=''){
  return(df_aux %>% knitr::kable(digits = ndec,caption = cap) %>% kableExtra::kable_styling())
}

Estudo

Novo relatório com análises para Ana Júlia - Los Hermanos

IRGA

UFSCar

CCSM

Rhizotron

dframe_rhiz = data.table::fread("/Users/rafaelcatoia/MyDrive/08_Consultorias/01_Ana_Julia/04_Projeto4/Dados/Rizhotron.txt")

dframe_rhiz_ufscar = dframe_rhiz %>%
  filter(Intituto=='UFSCar') %>% select(-Manejo,-Tratamentos) %>% 
  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')),
         Rhizotron_2022=10*Rhizotron_2022) %>% 
  pivot_longer(cols = c("Rhizotron_2021","Rhizotron_2022"),names_to = 'Safra',
               values_to = 'Rhizotron') %>% 
  mutate(Safra=gsub('Rhizotron_','',Safra))

dframe_rhiz_ccsm = dframe_rhiz %>%
  filter(Intituto=='CCSM') %>% select(-Tratamento,-Tratamentos) %>%  
  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV'))) %>% 
  pivot_longer(cols = c("Rhizotron_2021","Rhizotron_2022"),names_to = 'Safra',
               values_to = 'Rhizotron') %>% 
  mutate(Safra=gsub('Rhizotron_','',Safra))

UFSCar

dframe_summ_ufscar = dframe_rhiz_ufscar %>%
    select(Safra,RegHidrico,Tratamento,Rhizotron) %>%
    group_by(Safra,RegHidrico,Tratamento) %>% 
    summarise(Mean = mean(Rhizotron,na.rm=T),
              SD = sd(Rhizotron,na.rm=T)) %>% 
  mutate(`Mean ± SD`=paste(round(Mean,5),'±',round(SD,5)))

Graficos

Gráfico De Barras
Option 1
yname = 'Média de rhizotron'
dframe_summ_ufscar %>% ggplot(
      aes(x=Tratamento,
          fill=RegHidrico,
          y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    facet_wrap(~Safra)+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(yname)

Option 2
dframe_summ_ufscar %>% ggplot(
      aes(x=RegHidrico,
          fill=Tratamento,
          y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    facet_wrap(~Safra)+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(yname)

Gráfico De Barras com erro
Option 1
dframe_summ_ufscar %>% ggplot(
  aes(x=Tratamento,
      fill=RegHidrico,
      y=Mean)) +
  geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
  geom_errorbar(
    aes(x=Tratamento,y=Mean,
        color=RegHidrico,
        #linetype=RegHidrico,
        ymin=Mean-SD,ymax=Mean+SD),
    width = 0.2,alpha=0.8,stat = 'identity',color='black',
    position = position_dodge(0.9))+
  facet_wrap(~Safra)+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()+
  ylab(yname)

Option 2
dframe_summ_ufscar %>% ggplot(
  aes(x=RegHidrico,
      fill=Tratamento,
      y=Mean)) +
  geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
  geom_errorbar(
    aes(x=RegHidrico,y=Mean,
        color=Tratamento,
        #linetype=RegHidrico,
        ymin=Mean-SD,ymax=Mean+SD),
    width = 0.2,alpha=0.8,stat = 'identity',color='black',
    position = position_dodge(0.9))+
  facet_wrap(~Safra)+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()+
  ylab(yname)

Gráfico de Pontos + Média e Erro
Option 1
dframe_rhiz_ufscar %>% ggplot()+
    facet_wrap(~Safra,scales='free')+
    geom_point(aes(x=Tratamento,y=Rhizotron,
                          color=RegHidrico,
                          shape=RegHidrico),
               position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(x=Tratamento,y=Mean,
                          color=RegHidrico,
                          shape=RegHidrico),
               data=dframe_summ_ufscar,
               position = position_dodge(width = .75))+
    geom_errorbar(aes(x=Tratamento,y=Mean,
                          color=RegHidrico,
                          linetype=RegHidrico,
                      ymin=Mean-SD,ymax=Mean+SD),
                  data=dframe_summ_ufscar,
                  width = 0.2,alpha=0.9,
                  position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()

Option 2
dframe_rhiz_ufscar %>% ggplot()+
    facet_wrap(~Safra,scales='free')+
    geom_point(aes(x=RegHidrico,y=Rhizotron,
                          color=Tratamento,
                          shape=Tratamento),
               position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(x=RegHidrico,y=Mean,
                          color=Tratamento,
                          shape=Tratamento),
               data=dframe_summ_ufscar,
               position = position_dodge(width = .75))+
    geom_errorbar(aes(x=RegHidrico,y=Mean,
                          color=Tratamento,
                          linetype=Tratamento,
                      ymin=Mean-SD,ymax=Mean+SD),
                  data=dframe_summ_ufscar,
                  width = 0.2,alpha=0.9,
                  position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()

Tabela Resumo
dframe_summ_ufscar %>% select(Safra,RegHidrico,Tratamento,`Mean ± SD`) %>% kable_cust(cap = 'Média e Desvio Padrao - Rhizotron')
Média e Desvio Padrao - Rhizotron
Safra RegHidrico Tratamento Mean ± SD
2021 S BRF_PF+FD 0.01001 ± 0.00631
2021 S H10+FD 0.0105 ± 0.00623
2021 S H10+H152 0.03834 ± 0.01675
2021 TRF BRF_PF+FD 0.02045 ± 0.02714
2021 TRF H10+FD 0.00918 ± 0.01088
2021 TRF H10+H152 0.01859 ± 0.0209
2021 TRV BRF_PF+FD 0.00808 ± 0.00513
2021 TRV H10+FD 0.01402 ± 0.0136
2021 TRV H10+H152 0.02275 ± 0.00243
2022 S BRF_PF+FD 0.03457 ± 0.00308
2022 S H10+FD 0.03134 ± 0.00449
2022 S H10+H152 0.04866 ± 0.00287
2022 TRF BRF_PF+FD 0.02024 ± 0.00499
2022 TRF H10+FD 0.01622 ± 0.00192
2022 TRF H10+H152 0.02094 ± 0.00723
2022 TRV BRF_PF+FD 0.01795 ± 0.00295
2022 TRV H10+FD 0.02425 ± 0.00266
2022 TRV H10+H152 0.02961 ± 0.00519

Modelo

2021
df_tk <- dframe_rhiz_ufscar %>% filter(Safra=='2021')
fit_lm <- lm(data=df_tk,
          formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)

fit_aov <- aov(data=df_tk,
          formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)

summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)

summ_tk <- df_tk %>% group_by(RegHidrico,Tratamento) %>% 
  summarise(mean=mean(Rhizotron)) %>% arrange(-mean) 

summ_tk$Group =as.data.frame.list(mtcomp$`RegHidrico:Tratamento`)[,1]

summ_tk %>% kable_cust()
RegHidrico Tratamento mean Group
S H10+H152 0.03834 a
TRV H10+H152 0.02275 a
TRF BRF_PF+FD 0.02045 a
TRF H10+H152 0.01859 a
TRV H10+FD 0.01402 a
S H10+FD 0.01050 a
S BRF_PF+FD 0.01001 a
TRF H10+FD 0.00918 a
TRV BRF_PF+FD 0.00808 a
#tx <- with(data = df_tk,interaction(RegHidrico,Tratamento))
#amod <- aov(Rhizotron ~ RegHidrico+Tratamento+tx, data=df_tk)
#Tukey = agricolae::HSD.test(amod,trt = 'tx',group=TRUE)
#Tukey$groups %>% kable_cust()
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))
summ_aov
##                       Df   Sum Sq   Mean Sq F value Pr(>F)  
## RegHidrico             2 0.000107 0.0000534   0.257 0.7761  
## Tratamento             2 0.001277 0.0006387   3.077 0.0709 .
## RegHidrico:Tratamento  4 0.000847 0.0002116   1.020 0.4236  
## Residuals             18 0.003736 0.0002076                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.01001 0.00832 1.20307 0.24454
RegHidricoTRF 0.01044 0.01176 0.88738 0.38658
RegHidricoTRV -0.00192 0.01176 -0.16355 0.87191
TratamentoH10+FD 0.00049 0.01176 0.04153 0.96733
TratamentoH10+H152 0.02833 0.01176 2.40859 0.02695
RegHidricoTRF:TratamentoH10+FD -0.01176 0.01664 -0.70673 0.48878
RegHidricoTRV:TratamentoH10+FD 0.00544 0.01664 0.32723 0.74727
RegHidricoTRF:TratamentoH10+H152 -0.03019 0.01664 -1.81452 0.08630
RegHidricoTRV:TratamentoH10+H152 -0.01367 0.01664 -0.82149 0.42211
2022
df_tk <- dframe_rhiz_ufscar %>% filter(Safra=='2022')
fit_lm <- lm(data=df_tk,
          formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)

fit_aov <- aov(data=df_tk,
          formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)

summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)

summ_tk <- df_tk %>% group_by(RegHidrico,Tratamento) %>% 
  summarise(mean=mean(Rhizotron)) %>% arrange(-mean) 

summ_tk$Group = as.data.frame.list(mtcomp$`RegHidrico:Tratamento`)[,1]

summ_tk %>% kable_cust()
RegHidrico Tratamento mean Group
S H10+H152 0.04866 a
S BRF_PF+FD 0.03457 b
S H10+FD 0.03134 bc
TRV H10+H152 0.02961 bcd
TRV H10+FD 0.02425 bcde
TRF H10+H152 0.02094 cde
TRF BRF_PF+FD 0.02024 cde
TRV BRF_PF+FD 0.01795 de
TRF H10+FD 0.01622 e
#tx <- with(data = df_tk,interaction(RegHidrico,Tratamento))
#amod <- aov(Rhizotron ~ RegHidrico+Tratamento+tx, data=df_tk)
#
#Tukey = agricolae::HSD.test(amod,trt = 'tx',group=TRUE)
#Tukey$groups %>% kable_cust()
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))
summ_aov
##                       Df    Sum Sq   Mean Sq F value   Pr(>F)    
## RegHidrico             2 0.0017678 0.0008839  49.258 5.01e-08 ***
## Tratamento             2 0.0004839 0.0002419  13.483 0.000264 ***
## RegHidrico:Tratamento  4 0.0002684 0.0000671   3.739 0.021999 *  
## Residuals             18 0.0003230 0.0000179                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.03457 0.00245 14.13354 0.00000
RegHidricoTRF -0.01433 0.00346 -4.14327 0.00061
RegHidricoTRV -0.01661 0.00346 -4.80345 0.00014
TratamentoH10+FD -0.00323 0.00346 -0.93337 0.36298
TratamentoH10+H152 0.01409 0.00346 4.07497 0.00071
RegHidricoTRF:TratamentoH10+FD -0.00079 0.00489 -0.16097 0.87391
RegHidricoTRV:TratamentoH10+FD 0.00953 0.00489 1.94779 0.06721
RegHidricoTRF:TratamentoH10+H152 -0.01339 0.00489 -2.73656 0.01356
RegHidricoTRV:TratamentoH10+H152 -0.00244 0.00489 -0.49902 0.62381

CCSM

dframe_summ_ccsm = dframe_rhiz_ccsm %>%
    select(Safra,Manejo,Copa,PortaEnxerto,Rhizotron) %>%
    group_by(Safra,Manejo,Copa,PortaEnxerto) %>% 
    summarise(Mean = mean(Rhizotron,na.rm=T),
              SD = sd(Rhizotron,na.rm=T)) %>% 
  mutate(`Mean ± SD`=paste(round(Mean,5),'±',round(SD,5)))

Graficos

Gráfico de Barras
Option 1
yname = 'Média de rhizotron'
dframe_summ_ccsm %>% ggplot(
      aes(x=PortaEnxerto,
          fill=Copa,
          y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    facet_wrap(Safra~Manejo,nrow=1)+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(yname)

Option 2
dframe_summ_ccsm %>% ggplot(
      aes(x=Copa,
          fill=PortaEnxerto,
          y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    facet_wrap(Safra~Manejo,nrow=1)+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(yname)

Gráfico de Barras com erro
Option 1
dframe_summ_ccsm %>% ggplot(
  aes(x=PortaEnxerto,
          fill=Copa,
      y=Mean)) +
  geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
  geom_errorbar(
    aes(x=PortaEnxerto,y=Mean,
        color=Copa,
        ymin=Mean-SD,ymax=Mean+SD),
    width = 0.2,alpha=0.8,stat = 'identity',color='black',
    position = position_dodge(0.9))+
  facet_wrap(Safra~Manejo,nrow=1)+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()+
  ylab(yname)

Option 2
dframe_summ_ccsm %>% ggplot(
  aes(x=Copa,
          fill=PortaEnxerto,
      y=Mean)) +
  geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
  geom_errorbar(
    aes(x=Copa,y=Mean,
        color=PortaEnxerto,
        ymin=Mean-SD,ymax=Mean+SD),
    width = 0.2,alpha=0.8,stat = 'identity',color='black',
    position = position_dodge(0.9))+
  facet_wrap(Safra~Manejo,nrow=1)+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()+
  ylab(yname)

Gráfico de Pontos + Média e Erro
Option 1
dframe_rhiz_ccsm %>% ggplot()+
  facet_wrap(~Safra)+
  geom_point(aes(x=PortaEnxerto,
                 y=Rhizotron,
                 color=Copa,
                 shape=Manejo),
             position = position_dodge(width = .75),alpha=0.5)+
  geom_point(aes(x=PortaEnxerto,
                 y=Mean,
                 color=Copa,
                 shape=Manejo),
             data=dframe_summ_ccsm,
             position = position_dodge(width = .75))+
  geom_errorbar(aes(x=PortaEnxerto,
                 y=Mean,
                 color=Copa,
                 shape=Manejo,
                ymin=Mean-SD,ymax=Mean+SD),
                data=dframe_summ_ccsm,
                width = 0.2,alpha=0.9,
                position = position_dodge(0.75))+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()

Option 2
dframe_rhiz_ccsm %>% ggplot()+
  facet_wrap(~Safra)+
  geom_point(aes(x=Copa,
                 y=Rhizotron,
                 color=PortaEnxerto,
                 shape=Manejo),
             position = position_dodge(width = .75),alpha=0.5)+
  geom_point(aes(x=Copa,
                 y=Mean,
                 color=PortaEnxerto,
                 shape=Manejo),
             data=dframe_summ_ccsm,
             position = position_dodge(width = .75))+
  geom_errorbar(aes(x=Copa,
                 y=Mean,
                 color=PortaEnxerto,
                 shape=Manejo,
                ymin=Mean-SD,ymax=Mean+SD),
                data=dframe_summ_ccsm,
                width = 0.2,alpha=0.9,
                position = position_dodge(0.75))+
  theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
  scale_color_gdocs()+
  scale_fill_gdocs()

Tabela Resumo
dframe_summ_ccsm %>% select(Safra,Manejo,PortaEnxerto,Copa,`Mean ± SD`) %>% kable_cust(cap = 'Média e Desvio Padrao - Rhizotron')
Média e Desvio Padrao - Rhizotron
Safra Manejo PortaEnxerto Copa Mean ± SD
2021 Convencional FD BRF_PF 0.00521 ± 0.00401
2021 Convencional H152 BRF_PF 0.00837 ± 0.0044
2021 Convencional FD IAC10 0.0035 ± 0.00158
2021 Convencional H152 IAC10 0.01423 ± 0.00303
2021 Convencional FD IAC5 0.01404 ± 0.01049
2021 Convencional H152 IAC5 0.00987 ± 0.00401
2021 Ecologico FD BRF_PF 0.0088 ± 0.0014
2021 Ecologico H152 BRF_PF 0.01557 ± 0.00379
2021 Ecologico FD IAC10 0.00549 ± 0.00033
2021 Ecologico H152 IAC10 0.03573 ± 0.01341
2021 Ecologico FD IAC5 0.01022 ± 0.01135
2021 Ecologico H152 IAC5 0.01021 ± 0.00289
2022 Convencional FD BRF_PF 0.01011 ± 0.00569
2022 Convencional H152 BRF_PF 0.02145 ± 0.00038
2022 Convencional FD IAC10 0.00386 ± 0.00174
2022 Convencional H152 IAC10 0.00909 ± 0.00083
2022 Convencional FD IAC5 0.00171 ± 6e-05
2022 Convencional H152 IAC5 0.01492 ± 0.00074
2022 Ecologico FD BRF_PF 0.00425 ± 0.00071
2022 Ecologico H152 BRF_PF 0.02376 ± 0.00298
2022 Ecologico FD IAC10 0.00566 ± 0.00022
2022 Ecologico H152 IAC10 0.00891 ± 0.00142
2022 Ecologico FD IAC5 0.00459 ± 0.00013
2022 Ecologico H152 IAC5 0.01106 ± 0.00222

Modelo

2021
df_tk <- dframe_rhiz_ccsm %>% filter(Safra=='2021')
fit_lm <- lm(data=df_tk,
          formula = Rhizotron~Manejo+Copa +PortaEnxerto+
            Manejo*Copa+Manejo*PortaEnxerto +
            Copa*PortaEnxerto+
            Manejo*Copa*PortaEnxerto
          )

fit_aov <- aov(data=df_tk,
          formula = Rhizotron~Manejo+Copa +PortaEnxerto+
            Manejo*Copa+#Manejo*PortaEnxerto +
            Copa*PortaEnxerto+#+
            Manejo*Copa*PortaEnxerto
          )

summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)

summ_tk <- df_tk %>% group_by(Manejo,Copa,PortaEnxerto) %>% 
  summarise(mean=mean(Rhizotron)) %>% arrange(-mean) 

summ_tk$Group =as.data.frame.list(mtcomp$`Manejo:Copa:PortaEnxerto`)[,1]

summ_tk %>% kable_cust()
Manejo Copa PortaEnxerto mean Group
Ecologico IAC10 H152 0.03573 a
Ecologico BRF_PF H152 0.01557 b
Convencional IAC10 H152 0.01423 b
Convencional IAC5 FD 0.01404 b
Ecologico IAC5 FD 0.01022 b
Ecologico IAC5 H152 0.01021 b
Convencional IAC5 H152 0.00987 b
Ecologico BRF_PF FD 0.00880 b
Convencional BRF_PF H152 0.00837 b
Ecologico IAC10 FD 0.00549 b
Convencional BRF_PF FD 0.00521 b
Convencional IAC10 FD 0.00350 b
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))
summ_aov
##                          Df    Sum Sq   Mean Sq F value   Pr(>F)    
## Manejo                    1 0.0002372 0.0002372   5.620 0.026137 *  
## Copa                      2 0.0001736 0.0000868   2.056 0.149916    
## PortaEnxerto              1 0.0005457 0.0005457  12.928 0.001454 ** 
## Manejo:Copa               2 0.0002731 0.0001366   3.236 0.057000 .  
## Copa:PortaEnxerto         2 0.0007999 0.0004000   9.476 0.000926 ***
## Manejo:PortaEnxerto       1 0.0001861 0.0001861   4.409 0.046447 *  
## Manejo:Copa:PortaEnxerto  2 0.0001221 0.0000611   1.447 0.255160    
## Residuals                24 0.0010130 0.0000422                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.00521 0.00375 1.38946 0.17745
ManejoEcologico 0.00359 0.00530 0.67625 0.50535
CopaIAC10 -0.00171 0.00530 -0.32305 0.74946
CopaIAC5 0.00883 0.00530 1.66454 0.10901
PortaEnxertoH152 0.00316 0.00530 0.59530 0.55722
ManejoEcologico:CopaIAC10 -0.00160 0.00750 -0.21266 0.83339
ManejoEcologico:CopaIAC5 -0.00741 0.00750 -0.98735 0.33332
ManejoEcologico:PortaEnxertoH152 0.00362 0.00750 0.48201 0.63416
CopaIAC10:PortaEnxertoH152 0.00757 0.00750 1.00905 0.32302
CopaIAC5:PortaEnxertoH152 -0.00733 0.00750 -0.97650 0.33856
ManejoEcologico:CopaIAC10:PortaEnxertoH152 0.01589 0.01061 1.49799 0.14717
ManejoEcologico:CopaIAC5:PortaEnxertoH152 0.00054 0.01061 0.05116 0.95962
2022
df_tk <- dframe_rhiz_ccsm %>% filter(Safra=='2022')
fit_lm <- lm(data=df_tk,
          formula = Rhizotron~
            #Manejo+
            Copa +PortaEnxerto+
            #Manejo*Copa+Manejo*PortaEnxerto +
            Copa*PortaEnxerto#+
            #Manejo*Copa*PortaEnxerto
          )

fit_aov <- aov(data=df_tk,
          formula = Rhizotron~
            Manejo+Copa +PortaEnxerto+
            Manejo*Copa+Manejo*PortaEnxerto +
            Copa*PortaEnxerto+
            Manejo*Copa*PortaEnxerto
          )

summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)
tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)

summ_tk <- df_tk %>% group_by(Copa,PortaEnxerto) %>% 
  summarise(mean=mean(Rhizotron)) %>% arrange(-mean) 

summ_tk$Group =as.data.frame.list(mtcomp$`Copa:PortaEnxerto`)[,1]

summ_tk %>% kable_cust()
Copa PortaEnxerto mean Group
BRF_PF H152 0.02261 a
IAC5 H152 0.01299 b
IAC10 H152 0.00900 c
BRF_PF FD 0.00718 cd
IAC10 FD 0.00476 de
IAC5 FD 0.00315 e
par(mfrow=c(2,2))
plot(fit_lm)

par(mfrow=c(1,1))

par(mfrow=c(2,2))
plot(fit_aov)

par(mfrow=c(1,1))
summ_aov
##                          Df    Sum Sq   Mean Sq F value   Pr(>F)    
## Manejo                    1 0.0000021 0.0000021   0.480 0.495002    
## Copa                      2 0.0004488 0.0002244  50.641 2.44e-09 ***
## PortaEnxerto              1 0.0008701 0.0008701 196.364 4.73e-13 ***
## Manejo:Copa               2 0.0000099 0.0000050   1.122 0.342089    
## Manejo:PortaEnxerto       1 0.0000001 0.0000001   0.017 0.896199    
## Copa:PortaEnxerto         2 0.0001877 0.0000938  21.174 5.02e-06 ***
## Manejo:Copa:PortaEnxerto  2 0.0000871 0.0000435   9.827 0.000763 ***
## Residuals                24 0.0001063 0.0000044                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
summ_lm$coefficients %>% kable_cust()
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.00718 0.00107 6.71937 0.00000
CopaIAC10 -0.00242 0.00151 -1.60200 0.11963
CopaIAC5 -0.00403 0.00151 -2.66480 0.01228
PortaEnxertoH152 0.01543 0.00151 10.20594 0.00000
CopaIAC10:PortaEnxertoH152 -0.01118 0.00214 -5.23293 0.00001
CopaIAC5:PortaEnxertoH152 -0.00559 0.00214 -2.61646 0.01378

Potencial Hidrico

UFSCar

CCSM

Desenvolvimento Vegetativo

dframe_dv = data.table::fread("/Users/rafaelcatoia/MyDrive/08_Consultorias/01_Ana_Julia/04_Projeto4/Dados/DesenvVegetativo.txt")

#dframe_qf_ufscar = dframe_qf %>%
#  filter(Instituto=='UFSCar') %>% select(-Manejo) %>% 
#  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))
#
#dframe_qf_ccsm = dframe_qf %>%
#  filter(Instituto=='CCSM') %>% select(-Tratamento) %>%  
#  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))

UFSCar

CCSM

Qualidade do fruto

dframe_qf = data.table::fread("/Users/rafaelcatoia/MyDrive/08_Consultorias/01_Ana_Julia/04_Projeto4/Dados/QualidadeFruto.txt")

dframe_qf_ufscar = dframe_qf %>%
  filter(Instituto=='UFSCar') %>% select(-Manejo) %>% 
  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))

dframe_qf_ccsm = dframe_qf %>%
  filter(Instituto=='CCSM') %>% select(-Tratamento) %>%  
  mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))
summ_qf <- function(dframe,variable_name = ''){
  if (variable_name==''){
    stop()
  }
  dframe = dframe %>% select(one_of(variable_name),RegHidrico,Tratamento)
  
  
  colnames(dframe)[1]<-'Y'
  out<- list()
  
  dframe_summ <- dframe %>% group_by(RegHidrico,Tratamento) %>% 
    summarise(Mean = mean(Y,na.rm=T),
              SD = sd(Y,na.rm=T)) %>% 
    mutate(`Mean ± SD`=paste(round(Mean,2),'±',round(SD,2)))
  
  out$dframe_summ <- dframe_summ
  
  out$barp <- dframe_summ %>% ggplot(
    aes(x=RegHidrico,
        fill=Tratamento,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$barp_2 <- dframe_summ %>% ggplot(
    aes(x=Tratamento,
        fill=RegHidrico,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$barp_error <-  dframe_summ %>% ggplot(
    aes(x=RegHidrico,
        fill=Tratamento,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
    geom_errorbar(
      aes(x=RegHidrico,y=Mean,
          color=Tratamento,
          ymin=Mean-SD,ymax=Mean+SD),
      width = 0.2,alpha=0.8,stat = 'identity',color='black',
      position = position_dodge(0.9))+
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$barp_error_2 <-  dframe_summ %>% ggplot(
    aes(x=Tratamento,
        fill=RegHidrico,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
    geom_errorbar(
      aes(x=Tratamento,y=Mean,
          color=RegHidrico,
          ymin=Mean-SD,ymax=Mean+SD),
      width = 0.2,alpha=0.8,stat = 'identity',color='black',
      position = position_dodge(0.9))+
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$mean_desv_plot <-
    dframe %>% ggplot()+
    geom_point(aes(
      x=RegHidrico,
      y=Y,
      color=Tratamento,
      shape=Tratamento),
      position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(
      x=RegHidrico,
      y=Mean,
      color=Tratamento,
      shape=Tratamento),
      data=dframe_summ,
      position = position_dodge(width = .75))+
    geom_errorbar(aes(
      x=RegHidrico,
      y=Mean,
      color=Tratamento,
      shape=Tratamento,
      ymin=Mean-SD,ymax=Mean+SD),
      data=dframe_summ,
      width = 0.2,alpha=0.9,
      position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$mean_desv_plot_2 <-
    dframe %>% ggplot()+
    geom_point(aes(
      x=Tratamento,
      y=Y,
      color=RegHidrico,
      shape=RegHidrico),
      position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(
      x=Tratamento,
      y=Mean,
      color=RegHidrico,
      shape=RegHidrico),
      data=dframe_summ,
      position = position_dodge(width = .75))+
    geom_errorbar(aes(
      x=Tratamento,
      y=Mean,
      color=RegHidrico,
      shape=RegHidrico,
      ymin=Mean-SD,ymax=Mean+SD),
      data=dframe_summ,
      width = 0.2,alpha=0.9,
      position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))

  return(out)
}

summ_qf_ccsm <- function(dframe,variable_name = ''){
  if (variable_name==''){
    stop()
  }
  dframe = dframe %>% select(one_of(variable_name),Manejo,Copa,PortaEnxerto)
  
  colnames(dframe)[1]<-'Y'
  out<- list()
  
  dframe = dframe %>% mutate(Tratamento = paste(Copa,'+',PortaEnxerto,sep=''))
  
  dframe_summ <- dframe %>% group_by(Manejo,Tratamento) %>% 
    summarise(Mean = mean(Y,na.rm=T),
              SD = sd(Y,na.rm=T)) %>% 
    mutate(`Mean ± SD`=paste(round(Mean,2),'±',round(SD,2)))
  
  
  out$dframe_summ <- dframe_summ
  
  out$barp <- dframe_summ %>% ggplot(
    aes(x=Manejo,
        fill=Tratamento,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$barp_2 <- dframe_summ %>% ggplot(
    aes(x=Tratamento,
        fill=Manejo,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$barp_error <-  dframe_summ %>% ggplot(
    aes(x=Manejo,
        fill=Tratamento,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
    geom_errorbar(
      aes(x=Manejo,y=Mean,
          color=Tratamento,
          ymin=Mean-SD,ymax=Mean+SD),
      width = 0.2,alpha=0.8,stat = 'identity',color='black',
      position = position_dodge(0.9))+
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
   out$barp_error_2 <-  dframe_summ %>% ggplot(
    aes(x=Tratamento,
        fill=Manejo,
        y=Mean)) +
    geom_bar(stat = 'identity',position = 'dodge')+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
    geom_errorbar(
      aes(x=Tratamento,y=Mean,
          color=Manejo,
          ymin=Mean-SD,ymax=Mean+SD),
      width = 0.2,alpha=0.8,stat = 'identity',color='black',
      position = position_dodge(0.9))+
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$mean_desv_plot <-
    dframe %>% ggplot()+
    geom_point(aes(
      x=Manejo,
      y=Y,
      color=Tratamento,
      shape=Tratamento),
      position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(
      x=Manejo,
      y=Mean,
      color=Tratamento,
      shape=Tratamento),
      data=dframe_summ,
      position = position_dodge(width = .75))+
    geom_errorbar(aes(
      x=Manejo,
      y=Mean,
      color=Tratamento,
      shape=Tratamento,
      ymin=Mean-SD,ymax=Mean+SD),
      data=dframe_summ,
      width = 0.2,alpha=0.9,
      position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))
  
  out$mean_desv_plot_2 <-
    dframe %>% ggplot()+
    geom_point(aes(
      x=Tratamento,
      y=Y,
      color=Manejo,
      shape=Manejo),
      position = position_dodge(width = .75),alpha=0.5)+
    geom_point(aes(
      x=Tratamento,
      y=Mean,
      color=Manejo,
      shape=Manejo),
      data=dframe_summ,
      position = position_dodge(width = .75))+
    geom_errorbar(aes(
      x=Tratamento,
      y=Mean,
      color=Manejo,
      shape=Manejo,
      ymin=Mean-SD,ymax=Mean+SD),
      data=dframe_summ,
      width = 0.2,alpha=0.9,
      position = position_dodge(0.75))+
    theme_minimal(base_size=14)+theme(legend.position = 'bottom') + 
    scale_color_gdocs()+
    scale_fill_gdocs()+
    ylab(paste('Média de', variable_name ))

  return(out)
}

UFScar

Massa

results <- summ_qf(dframe_qf_ufscar,variable_name = "Massa")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 81.66667 7.09460 81.67 ± 7.09
S FD+H10 47.00000 2.64575 47 ± 2.65
S H14+H10 69.33333 3.51188 69.33 ± 3.51
S H152+H10 69.33333 6.11010 69.33 ± 6.11
S H228+H10 47.66667 0.57735 47.67 ± 0.58
TRF FD+BRS_PF 86.66667 8.14453 86.67 ± 8.14
TRF FD+H10 82.06667 4.10041 82.07 ± 4.1
TRF H14+H10 76.00000 2.64575 76 ± 2.65
TRF H152+H10 79.50000 8.50000 79.5 ± 8.5
TRF H228+H10 83.00000 8.18535 83 ± 8.19
TRV FD+BRS_PF 71.66667 1.15470 71.67 ± 1.15
TRV FD+H10 75.33333 1.15470 75.33 ± 1.15
TRV H14+H10 54.66667 8.50490 54.67 ± 8.5
TRV H152+H10 74.66667 3.51188 74.67 ± 3.51
TRV H228+H10 68.33333 7.02377 68.33 ± 7.02

Altura

results <- summ_qf(dframe_qf_ufscar,variable_name = "Altura")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 5.80000 0.20000 5.8 ± 0.2
S FD+H10 4.86667 0.11547 4.87 ± 0.12
S H14+H10 5.53333 0.05774 5.53 ± 0.06
S H152+H10 5.40000 0.20000 5.4 ± 0.2
S H228+H10 5.06667 0.05774 5.07 ± 0.06
TRF FD+BRS_PF 5.86667 0.28868 5.87 ± 0.29
TRF FD+H10 5.83333 0.15275 5.83 ± 0.15
TRF H14+H10 5.60000 0.00000 5.6 ± 0
TRF H152+H10 5.70000 0.10000 5.7 ± 0.1
TRF H228+H10 5.73333 0.25166 5.73 ± 0.25
TRV FD+BRS_PF 5.43333 0.11547 5.43 ± 0.12
TRV FD+H10 5.50000 0.10000 5.5 ± 0.1
TRV H14+H10 5.13333 0.11547 5.13 ± 0.12
TRV H152+H10 5.46667 0.32146 5.47 ± 0.32
TRV H228+H10 5.53333 0.28868 5.53 ± 0.29

Largura

results <- summ_qf(dframe_qf_ufscar,variable_name = "Largura")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 5.26667 0.11547 5.27 ± 0.12
S FD+H10 4.36667 0.05774 4.37 ± 0.06
S H14+H10 5.00000 0.10000 5 ± 0.1
S H152+H10 5.03333 0.15275 5.03 ± 0.15
S H228+H10 4.40000 0.10000 4.4 ± 0.1
TRF FD+BRS_PF 5.43333 0.15275 5.43 ± 0.15
TRF FD+H10 5.33333 0.11547 5.33 ± 0.12
TRF H14+H10 5.23333 0.05774 5.23 ± 0.06
TRF H152+H10 5.40000 0.10000 5.4 ± 0.1
TRF H228+H10 5.30000 0.17321 5.3 ± 0.17
TRV FD+BRS_PF 5.10000 0.00000 5.1 ± 0
TRV FD+H10 5.16667 0.05774 5.17 ± 0.06
TRV H14+H10 4.63333 0.30551 4.63 ± 0.31
TRV H152+H10 4.93333 0.28868 4.93 ± 0.29
TRV H228+H10 4.96667 0.15275 4.97 ± 0.15

AL

results <- summ_qf(dframe_qf_ufscar,variable_name = "AL")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 1.10333 0.02082 1.1 ± 0.02
S FD+H10 1.11667 0.02517 1.12 ± 0.03
S H14+H10 1.10667 0.01155 1.11 ± 0.01
S H152+H10 1.07333 0.03055 1.07 ± 0.03
S H228+H10 1.15333 0.04041 1.15 ± 0.04
TRF FD+BRS_PF 1.08333 0.02517 1.08 ± 0.03
TRF FD+H10 1.09333 0.02082 1.09 ± 0.02
TRF H14+H10 1.07333 0.01155 1.07 ± 0.01
TRF H152+H10 1.05667 0.00577 1.06 ± 0.01
TRF H228+H10 1.08333 0.02082 1.08 ± 0.02
TRV FD+BRS_PF 1.06667 0.02309 1.07 ± 0.02
TRV FD+H10 1.06667 0.02309 1.07 ± 0.02
TRV H14+H10 1.11000 0.05000 1.11 ± 0.05
TRV H152+H10 1.11000 0.01000 1.11 ± 0.01
TRV H228+H10 1.11333 0.03055 1.11 ± 0.03

Suco

results <- summ_qf(dframe_qf_ufscar,variable_name = "Suco")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 48.76667 4.05010 48.77 ± 4.05
S FD+H10 37.46667 3.35012 37.47 ± 3.35
S H14+H10 45.63333 1.50444 45.63 ± 1.5
S H152+H10 45.13333 2.58908 45.13 ± 2.59
S H228+H10 39.86667 0.46188 39.87 ± 0.46
TRF FD+BRS_PF 51.83333 2.10792 51.83 ± 2.11
TRF FD+H10 51.20000 2.07846 51.2 ± 2.08
TRF H14+H10 46.46667 1.76163 46.47 ± 1.76
TRF H152+H10 45.00000 2.70000 45 ± 2.7
TRF H228+H10 45.63333 3.13900 45.63 ± 3.14
TRV FD+BRS_PF 45.13333 0.05774 45.13 ± 0.06
TRV FD+H10 49.56667 0.75056 49.57 ± 0.75
TRV H14+H10 38.63333 4.06735 38.63 ± 4.07
TRV H152+H10 41.90000 1.60935 41.9 ± 1.61
TRV H228+H10 43.43333 3.76475 43.43 ± 3.76

Acidez

results <- summ_qf(dframe_qf_ufscar,variable_name = "Acidez")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 5.29067 0.38493 5.29 ± 0.38
S FD+H10 5.42190 0.28282 5.42 ± 0.28
S H14+H10 5.01870 0.39804 5.02 ± 0.4
S H152+H10 4.60677 0.51469 4.61 ± 0.51
S H228+H10 5.06987 0.18971 5.07 ± 0.19
TRF FD+BRS_PF 5.46163 0.26236 5.46 ± 0.26
TRF FD+H10 5.14260 0.43562 5.14 ± 0.44
TRF H14+H10 4.97363 0.08585 4.97 ± 0.09
TRF H152+H10 5.16717 0.15245 5.17 ± 0.15
TRF H228+H10 5.51613 0.34844 5.52 ± 0.35
TRV FD+BRS_PF 5.15400 0.32922 5.15 ± 0.33
TRV FD+H10 4.99920 0.20471 5 ± 0.2
TRV H14+H10 5.00657 0.06020 5.01 ± 0.06
TRV H152+H10 4.79327 0.15489 4.79 ± 0.15
TRV H228+H10 5.50537 0.33448 5.51 ± 0.33

BrixCorr

results <- summ_qf(dframe_qf_ufscar,variable_name = "BrixCorr")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 8.90000 0.17321 8.9 ± 0.17
S FD+H10 8.80000 0.10000 8.8 ± 0.1
S H14+H10 8.83333 0.15275 8.83 ± 0.15
S H152+H10 8.63333 0.05774 8.63 ± 0.06
S H228+H10 9.56667 0.20817 9.57 ± 0.21
TRF FD+BRS_PF 9.13333 0.15275 9.13 ± 0.15
TRF FD+H10 9.16667 0.05774 9.17 ± 0.06
TRF H14+H10 8.93333 0.11547 8.93 ± 0.12
TRF H152+H10 8.90000 0.10000 8.9 ± 0.1
TRF H228+H10 8.90000 0.00000 8.9 ± 0
TRV FD+BRS_PF 8.83333 0.20817 8.83 ± 0.21
TRV FD+H10 9.20000 0.10000 9.2 ± 0.1
TRV H14+H10 8.86667 0.20817 8.87 ± 0.21
TRV H152+H10 8.66667 0.05774 8.67 ± 0.06
TRV H228+H10 8.96667 0.25166 8.97 ± 0.25

Ratio

results <- summ_qf(dframe_qf_ufscar,variable_name = "Ratio")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 1.70000 0.10000 1.7 ± 0.1
S FD+H10 1.60000 0.10000 1.6 ± 0.1
S H14+H10 1.73333 0.11547 1.73 ± 0.12
S H152+H10 1.90000 0.20000 1.9 ± 0.2
S H228+H10 1.90000 0.10000 1.9 ± 0.1
TRF FD+BRS_PF 1.66667 0.05774 1.67 ± 0.06
TRF FD+H10 1.76667 0.15275 1.77 ± 0.15
TRF H14+H10 1.80000 0.00000 1.8 ± 0
TRF H152+H10 1.75000 0.05000 1.75 ± 0.05
TRF H228+H10 1.60000 0.10000 1.6 ± 0.1
TRV FD+BRS_PF 1.73333 0.15275 1.73 ± 0.15
TRV FD+H10 1.86667 0.05774 1.87 ± 0.06
TRV H14+H10 1.76667 0.05774 1.77 ± 0.06
TRV H152+H10 1.80000 0.10000 1.8 ± 0.1
TRV H228+H10 1.63333 0.05774 1.63 ± 0.06

IT

results <- summ_qf(dframe_qf_ufscar,variable_name = "IT")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 1.77000 0.18358 1.77 ± 0.18
S FD+H10 1.34333 0.11015 1.34 ± 0.11
S H14+H10 1.64667 0.03512 1.65 ± 0.04
S H152+H10 1.58667 0.09713 1.59 ± 0.1
S H228+H10 1.55333 0.02309 1.55 ± 0.02
TRF FD+BRS_PF 1.93000 0.05292 1.93 ± 0.05
TRF FD+H10 1.91667 0.05508 1.92 ± 0.06
TRF H14+H10 1.69000 0.08185 1.69 ± 0.08
TRF H152+H10 1.63667 0.08505 1.64 ± 0.09
TRF H228+H10 1.66000 0.11533 1.66 ± 0.12
TRV FD+BRS_PF 1.62667 0.04509 1.63 ± 0.05
TRV FD+H10 1.86000 0.02000 1.86 ± 0.02
TRV H14+H10 1.39667 0.18148 1.4 ± 0.18
TRV H152+H10 1.48000 0.06083 1.48 ± 0.06
TRV H228+H10 1.60000 0.17349 1.6 ± 0.17

VitC

results <- summ_qf(dframe_qf_ufscar,variable_name = "VitC")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
RegHidrico Tratamento Mean SD Mean ± SD
S FD+BRS_PF 28.54000 0.60622 28.54 ± 0.61
S FD+H10 35.68000 1.14582 35.68 ± 1.15
S H14+H10 32.45667 3.28649 32.46 ± 3.29
S H152+H10 37.78333 1.63102 37.78 ± 1.63
S H228+H10 33.36000 1.15763 33.36 ± 1.16
TRF FD+BRS_PF 33.21000 1.82371 33.21 ± 1.82
TRF FD+H10 32.76000 1.84984 32.76 ± 1.85
TRF H14+H10 32.50667 1.82870 32.51 ± 1.83
TRF H152+H10 34.24333 0.89131 34.24 ± 0.89
TRF H228+H10 34.38333 1.21187 34.38 ± 1.21
TRV FD+BRS_PF 33.45333 1.59080 33.45 ± 1.59
TRV FD+H10 35.30333 0.89891 35.3 ± 0.9
TRV H14+H10 38.17667 1.34500 38.18 ± 1.35
TRV H152+H10 37.89000 1.84789 37.89 ± 1.85
TRV H228+H10 35.40333 2.43097 35.4 ± 2.43

CCSM

Massa

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Massa")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 94.33333 11.50362 94.33 ± 11.5
Convencional BRF_PF+H152 89.33333 10.01665 89.33 ± 10.02
Convencional IAC10+FD 90.33333 4.72582 90.33 ± 4.73
Convencional IAC10+H152 86.66667 5.13160 86.67 ± 5.13
Convencional IAC5+FD 88.00000 12.12436 88 ± 12.12
Convencional IAC5+H152 81.66667 2.30940 81.67 ± 2.31
Ecologico BRF_PF+FD 77.33333 3.05505 77.33 ± 3.06
Ecologico BRF_PF+H152 83.66667 0.57735 83.67 ± 0.58
Ecologico IAC10+FD 95.00000 11.13553 95 ± 11.14
Ecologico IAC10+H152 90.33333 2.51661 90.33 ± 2.52
Ecologico IAC5+FD 80.33333 3.21455 80.33 ± 3.21
Ecologico IAC5+H152 88.00000 3.46410 88 ± 3.46

Altura

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Altura")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 6.10000 0.40000 6.1 ± 0.4
Convencional BRF_PF+H152 5.96667 0.25166 5.97 ± 0.25
Convencional IAC10+FD 6.06667 0.15275 6.07 ± 0.15
Convencional IAC10+H152 5.56667 0.11547 5.57 ± 0.12
Convencional IAC5+FD 5.90000 0.36056 5.9 ± 0.36
Convencional IAC5+H152 5.80000 0.10000 5.8 ± 0.1
Ecologico BRF_PF+FD 5.70000 0.10000 5.7 ± 0.1
Ecologico BRF_PF+H152 5.83333 0.05774 5.83 ± 0.06
Ecologico IAC10+FD 6.13333 0.15275 6.13 ± 0.15
Ecologico IAC10+H152 6.03333 0.11547 6.03 ± 0.12
Ecologico IAC5+FD 5.90000 0.10000 5.9 ± 0.1
Ecologico IAC5+H152 6.03333 0.05774 6.03 ± 0.06

Largura

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Largura")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 5.60000 0.36056 5.6 ± 0.36
Convencional BRF_PF+H152 5.40000 0.17321 5.4 ± 0.17
Convencional IAC10+FD 5.40000 0.10000 5.4 ± 0.1
Convencional IAC10+H152 5.10000 0.10000 5.1 ± 0.1
Convencional IAC5+FD 5.43333 0.25166 5.43 ± 0.25
Convencional IAC5+H152 5.26667 0.05774 5.27 ± 0.06
Ecologico BRF_PF+FD 5.16667 0.11547 5.17 ± 0.12
Ecologico BRF_PF+H152 5.26667 0.05774 5.27 ± 0.06
Ecologico IAC10+FD 5.56667 0.20817 5.57 ± 0.21
Ecologico IAC10+H152 5.23333 0.25166 5.23 ± 0.25
Ecologico IAC5+FD 5.30000 0.10000 5.3 ± 0.1
Ecologico IAC5+H152 5.43333 0.05774 5.43 ± 0.06

AL

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "AL")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 1.09000 0.01732 1.09 ± 0.02
Convencional BRF_PF+H152 1.10667 0.02082 1.11 ± 0.02
Convencional IAC10+FD 1.12333 0.01155 1.12 ± 0.01
Convencional IAC10+H152 1.09333 0.01155 1.09 ± 0.01
Convencional IAC5+FD 1.08667 0.02082 1.09 ± 0.02
Convencional IAC5+H152 1.10000 0.01000 1.1 ± 0.01
Ecologico BRF_PF+FD 1.10667 0.03055 1.11 ± 0.03
Ecologico BRF_PF+H152 1.10667 0.01528 1.11 ± 0.02
Ecologico IAC10+FD 1.10333 0.01155 1.1 ± 0.01
Ecologico IAC10+H152 1.15333 0.03786 1.15 ± 0.04
Ecologico IAC5+FD 1.11333 0.02082 1.11 ± 0.02
Ecologico IAC5+H152 1.11000 0.02000 1.11 ± 0.02

Suco

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Suco")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 50.00000 2.66646 50 ± 2.67
Convencional BRF_PF+H152 48.80000 1.57162 48.8 ± 1.57
Convencional IAC10+FD 48.30000 0.55678 48.3 ± 0.56
Convencional IAC10+H152 47.23333 2.61024 47.23 ± 2.61
Convencional IAC5+FD 48.76667 1.51438 48.77 ± 1.51
Convencional IAC5+H152 50.63333 1.25033 50.63 ± 1.25
Ecologico BRF_PF+FD 51.66667 2.72274 51.67 ± 2.72
Ecologico BRF_PF+H152 51.00000 0.34641 51 ± 0.35
Ecologico IAC10+FD 49.13333 1.68028 49.13 ± 1.68
Ecologico IAC10+H152 45.43333 3.35609 45.43 ± 3.36
Ecologico IAC5+FD 51.03333 1.51767 51.03 ± 1.52
Ecologico IAC5+H152 49.20000 1.76918 49.2 ± 1.77

Acidez

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Acidez")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 5.79277 0.21457 5.79 ± 0.21
Convencional BRF_PF+H152 4.73127 1.12089 4.73 ± 1.12
Convencional IAC10+FD 5.78470 0.14056 5.78 ± 0.14
Convencional IAC10+H152 4.64587 0.14255 4.65 ± 0.14
Convencional IAC5+FD 5.50670 0.55673 5.51 ± 0.56
Convencional IAC5+H152 5.28863 0.23447 5.29 ± 0.23
Ecologico BRF_PF+FD 5.43200 0.25991 5.43 ± 0.26
Ecologico BRF_PF+H152 5.71203 0.41669 5.71 ± 0.42
Ecologico IAC10+FD 5.43063 0.73447 5.43 ± 0.73
Ecologico IAC10+H152 5.57807 0.10602 5.58 ± 0.11
Ecologico IAC5+FD 5.59083 0.15735 5.59 ± 0.16
Ecologico IAC5+H152 5.64403 0.08699 5.64 ± 0.09

BrixCorr

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "BrixCorr")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plo 
## NULL
Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 9.03333 0.15275 9.03 ± 0.15
Convencional BRF_PF+H152 9.03333 0.15275 9.03 ± 0.15
Convencional IAC10+FD 9.06667 0.15275 9.07 ± 0.15
Convencional IAC10+H152 9.00000 0.30000 9 ± 0.3
Convencional IAC5+FD 9.26667 0.05774 9.27 ± 0.06
Convencional IAC5+H152 9.30000 0.17321 9.3 ± 0.17
Ecologico BRF_PF+FD 9.26667 0.15275 9.27 ± 0.15
Ecologico BRF_PF+H152 9.33333 0.15275 9.33 ± 0.15
Ecologico IAC10+FD 9.13333 0.15275 9.13 ± 0.15
Ecologico IAC10+H152 8.86667 0.05774 8.87 ± 0.06
Ecologico IAC5+FD 9.40000 0.10000 9.4 ± 0.1
Ecologico IAC5+H152 8.96667 0.15275 8.97 ± 0.15

Ratio

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Ratio")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 1.56667 0.05774 1.57 ± 0.06
Convencional BRF_PF+H152 1.96667 0.45092 1.97 ± 0.45
Convencional IAC10+FD 1.56667 0.05774 1.57 ± 0.06
Convencional IAC10+H152 1.93333 0.05774 1.93 ± 0.06
Convencional IAC5+FD 1.66667 0.15275 1.67 ± 0.15
Convencional IAC5+H152 1.76667 0.11547 1.77 ± 0.12
Ecologico BRF_PF+FD 1.73333 0.05774 1.73 ± 0.06
Ecologico BRF_PF+H152 1.63333 0.11547 1.63 ± 0.12
Ecologico IAC10+FD 1.70000 0.26458 1.7 ± 0.26
Ecologico IAC10+H152 1.56667 0.05774 1.57 ± 0.06
Ecologico IAC5+FD 1.66667 0.05774 1.67 ± 0.06
Ecologico IAC5+H152 1.56667 0.05774 1.57 ± 0.06

IT

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "IT")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 1.84667 0.06807 1.85 ± 0.07
Convencional BRF_PF+H152 1.80333 0.03215 1.8 ± 0.03
Convencional IAC10+FD 1.79000 0.02646 1.79 ± 0.03
Convencional IAC10+H152 1.73333 0.04726 1.73 ± 0.05
Convencional IAC5+FD 1.84000 0.05292 1.84 ± 0.05
Convencional IAC5+H152 1.92000 0.03606 1.92 ± 0.04
Ecologico BRF_PF+FD 1.94667 0.07506 1.95 ± 0.08
Ecologico BRF_PF+H152 1.94333 0.01528 1.94 ± 0.02
Ecologico IAC10+FD 1.83667 0.03215 1.84 ± 0.03
Ecologico IAC10+H152 1.64000 0.12124 1.64 ± 0.12
Ecologico IAC5+FD 1.95333 0.04163 1.95 ± 0.04
Ecologico IAC5+H152 1.79667 0.09018 1.8 ± 0.09

VitC

results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "VitC")
Gráfico de Barras
results$barp

results$barp_2

Gráfico De Barras com erro
results$barp_error

results$barp_error_2

Gráfico de Pontos + Média e Erro
results$mean_desv_plot

results$mean_desv_plot_2

Tabela Resumo
results$dframe_summ %>% kable_cust()
Manejo Tratamento Mean SD Mean ± SD
Convencional BRF_PF+FD 34.15667 2.19805 34.16 ± 2.2
Convencional BRF_PF+H152 37.52333 1.58418 37.52 ± 1.58
Convencional IAC10+FD 34.25333 1.33121 34.25 ± 1.33
Convencional IAC10+H152 33.71000 2.57874 33.71 ± 2.58
Convencional IAC5+FD 32.83667 2.87114 32.84 ± 2.87
Convencional IAC5+H152 33.98333 1.23678 33.98 ± 1.24
Ecologico BRF_PF+FD 31.96000 1.10585 31.96 ± 1.11
Ecologico BRF_PF+H152 33.32333 3.61777 33.32 ± 3.62
Ecologico IAC10+FD 34.23000 2.17690 34.23 ± 2.18
Ecologico IAC10+H152 36.66667 0.76846 36.67 ± 0.77
Ecologico IAC5+FD 31.53667 1.24340 31.54 ± 1.24
Ecologico IAC5+H152 33.23333 1.44901 33.23 ± 1.45